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An analytic method of smearing link variables in lattice QCD is proposed and tested. The differ- 
entiability of the smearing scheme with respect to the link variables permits the use of modern Monte 
Carlo updating methods based on molecular dynamics evolution for gauge-field actions constructed 
using such smeared links. In examining the smeared mean plaquette and the static quark-antiquark 
potential, no degradation in effectiveness is observed as compared to link smearing methods currently 
in use, although an increased sensitivity to the smearing parameter is found. 
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I. INTRODUCTION 

The extraction of hadron masses and matrix elements 
from Monte Carlo estimates of Euclidean-space corre- 
lation functions in lattice QCD can be done more reli- 
ably and accurately when operators which couple more 
strongly to the states of interest and less strongly to the 
higher-lying contaminating states are used. For states 
containing gluons, a crucial ingredient in constructing 
such operators in lattice QCD is link variable smear- 
ing. Operators constructed out of smeared or fuzzed 
links have dramatically reduced mixings with the high 
frequency modes of the theory. The use of such operators 
has been shown to especially benefit determinations of 
the glueball spectrumQ, hybrid meson masses0,0], the 
torelon spectrum^, and excitations of the static quark- 
antiquark potential|5|. 

Link variable smoothing is also playing an increasingly 
important role in the construction of improved lattice 
actions. In Ref. the use of so-called fat links in a 
staggered quark action was shown to significantly de- 
crease flavor symmetry breaking. Smeared links were 
subsequently used0 to construct hypercubic fermion ac- 
tions having improved rotational invariance. A staggered 
fermion action using a link smearing transformation's^ 
known as hypercubic (HYP) fat links was shown to im- 
prove flavor symmetry by an order of magnitude rela- 
tive to the standard action. The so-called Asqtad im- 
proved staggered quark action|^0,^ also makes use of 
link fattening to reduce flavor symmetry breaking. An- 
other variant of fermion actions which exploit smeared 
link variables is the fat link irrelevant clover (FLIC) 
action^^- FLIC fermions are Wilson-like and described 
by an action which includes an irrelevant clover improve- 
ment term constructed using smeared links. Fat links 
have also been used to construct a gauge action with 
reduced discretization errors using approximate rcnor- 
malization group transformations. 

The link fuzzing algorithm most often used in gluonic 
operator construction is that described in Ref. 14] in 
which every spatial link Uj{x) on the lattice is replaced 



by itself plus a real weight p times the sum of its four 
neighboring (spatial) staples, projected back into SU{2>). 
Such a fuzzing step is iterated rip times to obtain the 
final fuzzed link variables. Empirically, one finds that 
the projection into SU{i) is a crucial ingredient of the 
smearing. Link smearings which do not apply such a 
projection are found to be much less effective. The pro- 
jection into SU (3) is not unique and must be carefully 
defined so that all symmetry properties of the link vari- 
ables are preserved. Various ways of implementing this 
projection have appeared in the literature. Given a 3 x 3 
matrix V , its projection U into SU (3) is often taken to be 
the matrix U G SU{'i) which maximizes ReTr([/T/^'l'). An 
iterative procedure is required to perform such a maxi- 
mization. Alternatively the projected matrix U can 
be defined hyU {V^V)~^/^ Aei{V~^V'^Y/^ . 

Although the above smearing procedure works well in 
practice, the projection is somewhat unpalatable since it 
may be viewed as an arbitrary and abrupt way to remain 
within the group. More importantly, the branch cuts 
and lack of differentiability inherent in the projection can 
hinder or even make impossible the application of modern 
Monte Carlo updating techniques, such as hybrid Monte 
Carlo (HMC)[16|, which require knowing the response of 
the action to a small change in one of the link variables. 

A link smearing method which circumvents these prob- 
lems is proposed and tested in this paper. The link 
smearing method is analytic everywhere in the finite 
complex plane and utilizes the exponential function in 
such a manner to remain within SU{Z), eliminating the 
need for any projection back into the group. Because 
of this construction, the algorithm is useful for any Lie 
group. The method is described in Sec. Inland its prac- 
tical implementations when constructing operators and 
when computing the response of the action to a change 
in a link variable are detailed in Sees. IIIII and IIVI respec- 
tively. Numerical tests of the smearing are presented in 
Sec. Using the smeared plaquette and the effective 
mass aEcs{t) associated with the static quark-antiquark 
potential, no degradation in effectiveness is observed as 
compared to the standard link smearing method, al- 



though an increased sensitivity to the smearing param- 
eter is found. The resuhs also suggest that lattice ac- 
tions and operators constructed out of smeared link vari- 
ables may be much less afflicted by radiative corrections 
since the usually-dominant large tadpole contributions 
are drastically reduced. 



II. ANALYTIC LINK SMEARING 

A method of smearing link variables which is analytic, 
and hence differentiable, everywhere in the finite complex 
plane can be defined as follows. Let C^{x) denote the 
weighted sum of the perpendicular staples which begin 
at lattice site x and terminate at neighboring site x+fi: 

C^,{x) ^ ^p^^f [/^(a;)t/^(x-}-z))L/J;(a; + A) 

+Ul{x-0)Uf,{x-iy)U^{x-i)+fi)^, (1) 

where /t, u are vectors in directions n, v, respectively, hav- 
ing the length of one lattice spacing in that direction 
(fundamental lattice translation vectors). The weights 
p^i/ are tunable real parameters. Then the matrix Q^{x), 
defined in SU{N) by 

Q^{x) = '-in\Xx)-^l^{x^^^Tr{nl{x)~n^{x^, 
^fj,{x) = Cfj_{x) Uj^{x), (no summation over fi) (2) 

is Hermitian and traceless, and hence, 6**3 is an el- 
ement of SU{N). We use this fact to define an itera- 
tive, analytic link smearing algorithm in which the links 
ujr\x) at step n are mapped into links ulr^^\x) using 



C/("+i)(x)=exp(zg(,")(x)) U^^Ha 



(3) 



The fact that e^'^^Ca^) jg an element of SU{N) guarantees 

that ujr^^\x) is also an element of SU{N), eliminating 
the need for a projection back onto the gauge group. This 
fuzzing step can be iterated rip times to finally produce 
link variables which we call stout links\ll\. denoted by 



u. 



One common choice of the staple weights is 



(4) 



(5) 



which yields a three-dimensional scheme in which only 
the spatial links are smeared. Another common choice is 
an isotropic four-dimensional scheme in which all weights 
are chosen to be the sam e Pu i, — p. Note that such a 
scheme was used in Ref. |lg as a field transformation 
to eliminate a certain interaction term from the most 
general 0{a^) on-shell improved gauge action. 



FIG. 1: The expansion up to first order in tlie p^^ of the new 
link variable U^^^ in terms of paths of the original links. Each 
closed loop includes a trace with a factor 1/A'' in SU{N). 



It is not difhcult to show that, given an appropri- 
ate choice of weights p^^, the stout links have symme- 
try transformation properties identical to those of the 
original link variables. Under any local gauge transfor- 
mation G{x), the link variables transform according to 
Uf_i{x) G{x) Ufj.{x) G'^{x + [i), where the G{x) are 
SU{N) matrices. Thus, under such a gauge transfor- 
mation. 



G{x) Q^{x) G\x), 



G{x) 



4^} 



GHx), 



(6) 
(7) 



C///)(:e) ^ Gix) e'^^^^^G^ix)Gix)Uf,{x) G^{x + fi) 

= Gix) Ul'Hx) GHx+fi), (8) 

as required. One can also easily show that if the weights 
Pfif respect the rotation and reflection symmetries of the 
lattice, the stout links obey the required transformation 
properties under all rotations and all reflections in a plane 
containing the link. Lastly, consider a reflection in a 
plane normal to a link U^{x) and passing through its 
midpoint x-\-^ji. Under such an operation, the link Uf^{x) 
transforms according to Uf^{x) U-^{x+jj,) = Uj^{x) and 
for ^ /I, U„{x) ^ U,j{x + fi) and U},{x + jl) ^ U},{x), 
yielding G^{x) Cjj^{x) assuming the p^i, are real, so 
that 

n^{x) ^ Gl{x)U,,{x) = Ul{x)nl{x)U^{x), (9) 
Q^.{x) ^ -Ul{x)Q^{x)U,,{x), (10) 



U'-'Kx) ^ [/t(x)e-'«-(-)=C/Wt(^), 



(11) 



as required. Given that U^p {x) transforms under all 
symmetry operations in the same manner as C/^(x), it 
follows that the stout links Ufj_{x) have symmetry trans- 
formation properties identical to those of the original link 
variables. 

Since the exponential function has a power series ex- 
pansion with an infinite radius of convergence, each stout 
link may be viewed as an incredibly large and compli- 
cated sum of paths. For small p^i,, the paths which 
make up the link variable U^^^ to first order in the p^j^ 
are shown in Fig. ^ Note that the standard smearing 
method, defined by 



U^Pix)^Vsui3) \u,ix)+G^{x) 



(12) 
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where Vsu{3) denotes the projection mto SU{3), yields 
the same sum of paths at first order in p^,^. 

A few further remarks are worthy of note. First, an 
alternative smearing scheme in which the weights p^^ are 
chosen to be imaginary does not reduce the couplings to 
the high-frequency modes of the theory. Secondly, it is 
not possible to remove the leading 0{a^) discretization 
errors in the Wilson gauge action by expressing the action 
in terms of stout links; diminishing the 0{asa^) errors by 
reducing their tadpole contributions is the best which can 
be achieved. Thirdly, the use of stout links in Symanzik- 
improved actions may eliminate the need to use tadpole- 
improvement in determining the couplings in the action. 
Lastly, the definition of C^{x) given in Eq. is not 
unique. Analytic smearing schemes can be constructed 
as outlined above using sums of other paths from site x 
to x + jx, such as those used in fat hnk|6j and HYP,^] 
smearings. 



III. IMPLEMENTATION OF THE SMEARING 



numbers satisfying 91 + 92 + 93 = 0. These eigenvalues 
are the three roots of a cubic polynomial which can be 
easily determined: 



91 = 2u, 

92 = -u + w, 

93 = -91 - 92 = -u 



where 



u 
w 



= arccos 



|ci cos (i6') , 
ci sm ( i^) , 
Co 







(20) 
(21) 
(22) 



(23) 
(24) 

(25) 



Given these eigenvalues, the matrix Q can be written 



Q = M Aq M- 



Ao = 



91 
92 
(73 



(26) 



To implement this analytic link smearing scheme, the 
efficient evaluation of exp{iQ), where Q is a traceless Her- 
mitian 3x3 matrix, is required. The Cayley-Hamilton 
theorem states that every matrix is a zero of its charac- 
teristic polynomial, so that 



where 



Co = detQ= iTr(Q3), 
ci = iTr(g2)>0. 



(13) 

(14) 
(15) 



The Hermiticity of Q requires 27cq — 4cf < and the def- 
inition of Q given in Eq. restricts the possible values 
of ci. Thus, the coefficients cq and ci satisfy 



max 
Co 



where 



< Co < cS^^ 



< ci < c'l" 



max 
Cl 



1 

32 



(69 + ll\/33)(^p^. 



(16) 

(17) 
(18) 



for each /i. Eq. (|13() implies that Q" for integer n > 3 can 
be expressed in terms of Q^, Q, and the identity matrix 
/. Hence, we can write 



/o / + /i Q + /2 Q\ 



(19) 



where the three scalar coefficients fj = fj (co, ci) are basis 
independent, depending only on cq and ci. Eq. H19|l is 
valid for any 3x3 Hermitian, traceless matrix Q. 

Let (71,(72,93 denote the three eigenvalues of Q. Since 
Q is Hermitian and traceless, we know that these are real 



1 91 9? 










1 92 qj 




h 






1 93 93 




.f2 




e»93 



where Af is a unitary matrix. Then it easily follows that 

e^^« =/o/ + /iAq + /2A|. (27) 

Explicitly, we have the following linear system of equa- 
tions to solve: 



(28) 



If all three eigenvalues are distinct, this system of equa- 
tions has a unique solution, but when two of the eigen- 
values are exactly the same, the solution to Eq. (|28|l is 
not unique and one of the /j's can be freely set to any 
value. The case of degenerate eigenvalues occurs when 
Co ^ ±c^-^ 

In practice, it is extremely unlikely that an exact de- 
generacy will be encountered during a numerical simula- 
tion. Much more likely is the possibility that two eigen- 
values are nearly, but not quite, equal. Encounters with 
both near and exact degeneracies can be handled by ex- 
pressing the fj in terms of u and w in order to isolate 
and tame the numerically sensitive part. One finds that 
the fj factors can be written 



hi 



(91*2-^2)' 

where the hj are well-behaved functions given by 
ho = (M2-u,2)e2» -)-e-™{8w2cos(u;) 



(29) 



(30) 



hi = 2Me2™ - e~™{2ucos(u;) 

-^(3u2-^i;2)eo(^«)}, (31) 
h2 = e2*"~e-™{cos(w)-K3i<o(ii')}, (32) 
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defining 



C / \ ^^^^ /QQ^ 

CoH = . (33) 



w 



The w ^ problem as cq c™^^ has been completely 
contained in ^o- To numerically evaluate this factor, one 
uses, for example. 



sin(w)/w. 



{l-^w^)), |u;|<0.05, 
\w\ > 0.05. 



However, the w ^ 3u ^ v3/2 limit as co —c^'^^ has 
not been tamed. Fortunately, this problem can be cir- 
cumvented using the following symmetry relations under 
Co -Cq: 



/j(-co,ci) = i-y /*(co,ci). 



(34) 



Thus, the determination of the fj coefficients for cq < 
can be achieved by computing the coefficients for |co| and 
utiHzing Eq. This means that a situation in which 

the denominator in Eq. (|29|l becomes nearly zero will 
never be encountered. It should be emphasized that the 
fj coefficients are smooth, non-singular functions of cq 
and ci. There are no actual singularities as cq — > zLc™'^^. 
The numerical evaluation of the fj coefficients to machine 
precision in these limits simply requires additional care. 

As an aside, if one writes Q — ^ J2l=i Qa^'^ in terms 
of the eight Gell-Mann matrices A", then 





= '"0 + iELl"a''^°' 


(35) 


Uo 


= ./O + |ci/2, 


(36) 


Ua 




(37) 


Co 


"12^ Q aQbQ Cl 


(38) 


Cl 


— "^QaQai 


(39) 



where d and f are the real symmetric and antisym- 
metric structure constants, respectively. 



IV. MOLECULAR DYNAMICS EVOLUTION 

The analyticity of the stout-link scheme permits the 
use of modern Monte Carlo updating methods for gauge- 
field actions constructed using stout links. Molecular dy- 
namics evolution forms the core of many of the Markov 
processes used to generate ensembles of field coiifig^ura- 
tions needed for unquenched QCD simulations 16, 19, 20, 
|2]| . Since a small change to the underlying gauge fields 
always leads to a small change in the stout links, then 
provided the small change to the stout links causes small 
changes to the action of the theory, the force term on 
the underlying links is always well-defined. The impor- 
tance of defining a link smearing scheme which permits 
the computation of the molecular dynamics force term 
has been emphasized recently in Ref. p^. 



The S^(2:) force field describes how the action S 
changes when the gauge-field link variables U^{x) change, 
holding the momenta Pfi{x) and the pseudofermion field 
4>{x) fixed. For this reason, we shall write the action 
S = S[U] in this section, ignoring its dependence on the 
pseudofermion field. The transpose of the force field is 
defined by the derivative of the action with respect to the 
link variables: 



/ dS[U] 



(40) 



Assume that the action can be written as a sum of a 
term Sth [U] constructed entirely of the original thin link 
variables and a term Sst[U] constructed completely out 
of stout links: 

S[U]^Stb[U] + Sst[U]. (41) 
Then the force field may be written 

S^(x) = I](*h)(^)^j.(o)(^)^ (42) 

where 



The computation of S^*^^ (a;) is usually straightforward, 
and nothing more about its determination needs discus- 
sion here. We now focus on the evaluation of E^°''(x). 

Recall that the stout links are constructed iteratively 
starting with the original links U — > U'^^^ J7^^^ 
■ ■ ■ ~f [/("p) — U. The computation of the force field 
proceeds similarly in an iterative fashion, except that the 

order is reversed E = E^"") -> E'"""!) ^ > E^^) 

E^"^ The sequence starts by computing 



E^(x) 



(44) 



This step depends on the form of the gauge and fermion 
action in terms of the stout links and is usually just as 

straightforward as the computation of E^**^-* (x). In subse- 
quent steps, E^*"') and C/^^"^) are used to compute E^'^"^^ 
where the effective force at level k is defined by 



The recursive mapping 

|j;(/c)^^(fc-l)| ^ J](fe-1) 



(45) 



(46) 



is repeated until E^*^-* is finally evaluated. 

In Monte Carlo methods based on molecular dynamics, 
such as HMCllll and the i?-algorithm0, the flow of 



lattice configurations through phase space via Hamilton's 
equations is parametrized by a fictitious simulation time 
coordinate r. To determine the mapping in Eq. 14fci|l . it is 
convenient to express the force field in terms of the rate of 
change of the action with respect to this time coordinate 
r: 

|-5,t[C/] =25]ReTr|l]W(x) ^Ul'^H^)] , (47) 

for /c = 0, 1, . . . , Up. One then proceeds using the chain 
rule of differentiation: 



dr 



iQ 



dr 



d{e 



dr 



(48) 



To simplify matters, Q)!'~''{x), Ul:~''[x), E}; {x) 
shall simply be written asQ,U,Yi, respectively, in the cal- 
culations which follow, and y!"^^ (x) and ult^ (x) shall be 
written as E' and [/', respectively. The Cay ley-Hamilton 
theorem gives us 



d{e 



dr 



^ (/o + AO + /2O') , 
dr dr dr 



dO dO dO , , 

Since the fj coefficients are functions fj — fj{u, w) of u 
and w only, one has 



^^/df^\du^/df^\dM^^ 
dr \ du J dr \ dw J dr ' 



(50) 



and since u and w are functions of cq and ci only, then 

(51) 
(52) 



du 
dw 

Next, one finds that 



du \ dcQ / du \ dci 

dco J dr \ dci J dr ' 

dw \ dcQ / dw \ dci 

dcQ J dr \ dci J dr 



dco 



dQ 



,dQ 



(53) 



f ^lATr(g^)^Tr^g^), (54) 



and 



du 



du 



dcQ 2(9^2 — ■u;2) ' (9u2 — w"^) ' 



dw 



-3u 



dw 



(55) 



dcQ 2w(9m^ — w"^) ' dci 2w{9u'^ — w^) 



Now define 



,(1) ^ ^(2) ^ 

c)u ' w dw ' 



and 



2ur[^^ + (3^2-^2)^(2) _ 2(15^2+^2)/^. 
bi] = ^ 



2(9w2 _ u;2)2 



then 



(56) 



,(57) 



(58) 



5yTr ( Q-^ ] + 62, Tr [Q^j, (59) 



^ - Surf ^ - 24m/,- 



"2j 0/'n„,2 „„2^2 



2(9^2 -u;2)2 



_^5i,Tr 0-^ 



where 



with 



= 2(7. + i(^i2_u;2))e2» 

-H2e~™|4M(2-m) cosu; 

-Hi (97/2 + u)2 - m(3u2 + w2))^o(u;)|, (60) 
rj^' = 2(l + 2m)e2'« + e-*"|-2(l-m)cosu; 

+i{Qu + i{w'^ -■iu^))io{w)Y (61) 
r^^^ = 2ie2™-Hie-™|cosw-3(l-m)Co(w)},(62) 
r(2) ^ -2e2"'-F2me-™|cosu; 

+ (l-h4m)eo(«^) +37/2^1 (u)) I , (63) 
rf^ = -ie"™|cosw (1 2w)^o(w) 

'ei(^")}, (64) 

(65) 

(66) 
(67) 



J2) 



e-™ { Co (u;)-3z<i («;)}, 



6(w) = 



smui 

1 

w 

cos w sin w 



Given the above results. 



d(e^Q) 



= Trlg:^ )i3i + Tr(g2^\B2 



,dQ 
dr 



dr 



^,dQ dQ dQ 

+fi-^ + f2-^Q + f2Q-^, (68) 



where 



B, = 6,0 + ba Q + bi2 Q^ 



(69) 



The numerically sensitive w — s- limit has been totally 
absorbed into ^0 and ^1 , and numerical problems as co — *■ 
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can be circumvented by exploiting the following 



symmetry relation: 



6»,(-co,ci)-(-)'+^'+i b*Jco,c,). 



(70) 



Once again, we emphasize that these coefficients are well- 
behaved, non-singular functions of cq and ci. 

The rate of change of the stout links with respect to 
the simulation time is then given by 



dU' 



and hence, 



ar dr 



Bi + Tr [Q 



,dQ 
dr 



Bo 



(71) 



ReTr E 



^,dU^ 
dr 



^ReTrlEV'? — 
dr 



ReTr 



defining 



(72) 



A 



Vr + rt)-4Tr(r + rt), 



2 ' ' 2N' 

r = Tr(E'BiC/) Q + Tr(S'B2C/) 



(73) 



(74) 



One sees that t')^~^\x), and hence, h}i^~^'{x), are de- 
fined on each lattice link in terms of the link variables 
uIj' ^\x) and y!"^\x). At this point, the details of the 
staple construction in C^(a;) must be included. Given 
Eq. JQ) and utilizing trace cyclicity and translational in- 
variance, one eventually obtains the following recursion 
relation: 

i:^{x) - S;,(x)(/o/ + hQ^{x) + f2Ql{x)) 
+iCl{x)K^{x) - iY,[ 

Put,U^ix+fi)Ul{x + v)Ul{x)A^{x) 
+PiiiyUl{x~v+jl)Uj^{x-i')A^j_{x~v)Uiy{x-v) 
+PvtLUl{x-v+jl)K^{x-v+jl)U'l^{x-v)Uy{x-v) 
-PvtiUl{x~v+jji)Ul_{x~i')K^{x-v)Ui,{x-v) 
-p^^K{x + fi)U,,{x + (i)Ul{x+v)Ul{x) 

+p^J,^U^{x+p)Ul{x + v)^^,{x+v)Ul{x)Y (75) 

in which miprimed quantities refer to step (fc — 1) and 
primed quantities refer to step k. 

Although the force field is related to the underlying 
link variables in a complicated manner, each step in the 
recursive computation of the force field outlined above is 
straightforward, facilitating a natural and efficient imple- 
mentation in software. Analyticity in the entire (finite) 
complex plane is an important property here; if S^'^-' (x) 

is well behaved, then S^"-* (x) must be also. Note that the 
above formalism can be easily adapted for other defini- 
tions of C^{x). 



(fe-i) , 




FIG. 2: The efltective energy aE^s{r) defined in Eq. ^ 
for a static quark-antiquark pair separated by a distance 
r = 5a for several levels of smearing Up = 1,5, and 20 
against the smearing parameter p. These results were ob- 
tained on a 12* lattice using the Wilson gauge action with cou- 
pling P — 5.7. Curves labeled A„^ indicate results using the 
spatially-isotropic three-dimensional version of the analytic 
stout hnk smearing scheme with Up — 1,5, and 20 levels, while 
the curves labeled Bnp show the results for links smeared us- 
ing Eq. 11211 with the projection method of Ref. fisf . 



NUMERICAL TESTS 



The efficacy of stout links both as a gluonic-operator 
smearing algorithm and for reducing the effects of ultravi- 
olet gluon modes in short-distance quantities was tested 
in several Monte Carlo simulations. 

The energy of gluons in the presence of a static quark- 
antiquark pair separated by a distance r can be extracted 
from r X t Wilson loops W{r, t) as the temporal extent t 
becomes large. The Wilson loop can be viewed as the cor- 
relation function of a gauge-invariant operator consisting 
of a static quark-antiquark pair connected by a product 
of link variables following a straight-line path between 
the quark and the antiquark. Extraction of the lowest 
energy can be done much more reliably if the couplings 
of the quark-antiquark-gluon operator with higher-lying 
states are small. This can be achieved by utilizing a 
product of smeared links connecting the quark and the 
antiquark, instead of the original link variables. In other 
words, the ground state energy of gluons in the presence 
of a static quark-antiquark pair can be more reliably de- 
termined from Wilson loops W{r,t) constructed using 
smeared spatial links in the r-link paths on the initial 
and final time slices. A measure of how well the cou- 
plings to the higher levels are reduced is the effective 
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FIG. 3: The effective energy aEcs{r) defined in Eq. l|76)l 
for a static quark-antiquark pair separated by a distance 
r — 10a for several levels of smearing rip = 1,5, and 20 
against the smearing parameter p. Results were obtained 
on a 24* lattice using the Wilson gauge action with cou- 
pling /3 = 6.2. Curves labeled An^ indicate results using the 
spatially-isotropic three-dimensional version of the analytic 
stout link smearing scheme with Up = 1, 5, and 20 levels, while 
the curves labeled Bn^ show the results for links smeared us- 
ing Eq. II12II with the projection method of Ref. (Tsf. 

energy for a single time step, defined by 

Under unexceptional circumstances for a lattice gauge ac- 
tion with a positive-definite transfer matrix, such as the 
Wilson action, reducing the excited-state contamination 
tends to lower the effective energy given above. 

The effective energy defined in Eq. H76() was used to test 
the ability of the stout link smearing scheme to reduce 
excited-state contamination in gluonic operators. Results 
for r = 5a on a 12^ lattice using the Wilson gauge action 
with coupling /3 = 5.7 are shown in Fig. [2 and results for 
r = 10a on a larger 24"* lattice with /3 = 6.2 are shown in 
Fig. 121 The results are compared to those obtained with 
links smeared using Eq. H12|) and the projection method 
of Ref. [13. 

Both figures show that link smearing can dramatically 
reduce contamination from the high-lying modes of the 
theory. Consider first the results obtained using the 
spatially-isotropic three-dimensional version of the stout 
link smearing scheme {pjk = p, = = 0). For 
a given number of smearing iterations Up, the effective 
energy decreases initially as the smearing parameter p is 
increased from zero. Eventually an optimal value for the 
smearing parameter is reached at which the effective en- 
ergy is minimized. This optimal value decreases as the 




FIG. 4: The mean smeared plaquette against the smearing 
parameter p. These results were obtained using the Wilson 
gauge action on a 12* lattice with coupling /3 = 5.7. Curves 
labeled A„p indicate results obtained using the isotropic four- 
dimensional version of the analytic stout link smearing scheme 
with Up — 1 and 5 levels, while the curves labeled B„p show 
the results with links smeared using Eq. (1121 and the projec- 
tion method of Ref. [Tsf . 

number of smearing levels rip is increased, and the reduc- 
tion in the effective energy at this optimal p value is sub- 
stantial. Further increasing p beyond its optimal value 
then results in a sharply increasing effective energy. The 
rapidity of both the fall and rise of the effective energy 
about its minimum is more pronounced for larger Up. 
Also, the minimum value decreases as Up increases until 
eventually a saturation point is reached beyond which no 
additional reduction occurs. 

The results obtained with links smeared using Eq. l(T^ 
and the projection method of Ref. [l3| display essentially 
the same trends, except that each minimum in the effec- 
tive energy is much broader. The increased sensitivity 
of the stout link smearing scheme to the parameter p is 
not surprising since p occurs inside an exponential func- 
tion. However, it is important to note that both smear- 
ing methods produce nearly the same minimum values of 
the effective energy. The stout link smearing scheme is 
just as effective at reducing excited-state contamination 
in gluonic operators as the standard link fuzzing scheme 
in current use, although more careful tuning is necessary 
due to the increased sensitivity to the smearing param- 
eter p. Note that this analytic link smearing method 
has already been successfully applied in computing the 
spectrum of torelon excitations'4j. 

The main purpose of using smeared links in lattice 
gauge and fermion actions is to reduce discretization ef- 
fects caused by the ultraviolet modes in the lattice theory. 
Often such effects are dominated by large contributions 
from tadpole diagrams. A simple measure of how well a 
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link smearing scheme can reduce artifacts from tadpole 
diagrams is the mean smeared plaquette. A value of the 
mean smeared plaquette near unity indicates a substan- 
tial reduction of the tadpole contributions. 

Results for the mean smeared plaquette on a 12-* lattice 
using the Wilson gauge action with coupling /3 = 5.7 are 
shown in Fig. ^ For isotropic four-dimensional versions 
(P/if — p) of both smearing methods, the mean smeared 
plaquette initially increases towards unity as p increases 
from zero until a maximum is reached. The two smear- 
ing methods produce nearly the same maximum value. 
As p is increased further, the mean smeared plaquette in 
the stout link scheme quickly begins to fall, whereas the 
standard smearing scheme changes little. As np initially 
increases from zero, the maximum value of the mean pla- 
quette also increases, but eventually a saturation point 
is reached. Note that the maximum values for Up — b 
are very near to unity, suggesting a dramatic reduction 
of tadpole contributions. In summary, the analytic stout 
link smearing scheme is observed to be just as efficacious 
as the standard smearing scheme in reducing discretiza- 
tion effects from the ultraviolet gluon modes, but with 
an increased sensitivity to the smearing parameter p. 

VI. CONCLUSION 

Link- variable smoothing is a crucial ingredient in con- 
structing gluonic operators which have dramatically re- 
duced mixings with the high frequency modes of the the- 
ory. Link smearing is also playing an increasingly impor- 
tant role in the construction of improved lattice actions. 
The lack of differentiability with respect to the underly- 



ing link variables of standard smearing schemes prevents 
the use of efficient Monte Carlo updating methods based 
on molecular dynamics evolution. 

A link smearing method which circumvents these prob- 
lems was proposed and tested in this paper. The link 
smearing method is analytic everywhere in the finite com- 
plex plane and utilizes the exponential function in such a 
manner to remain within SU (3), eliminating the need for 
any projection back into the group. Because of this con- 
struction, the algorithm is also useful for any Lie group. 
An efficient implementation of this smearing scheme, as 
well as the recursive computation of the force term de- 
scribing the change of the action in response to a varia- 
tion of the link variables, was described. Although the 
force field is related to the underlying link variables in 
a very complicated manner, each step in the recursive 
computation of the force field is actually straightforward, 
allowing a natural and efficient design in software. The 
smeared mean plaquette and the effective energy associ- 
ated with the static quark-antiquark potential were used 
to show that no degradation in effectiveness is observed 
as compared to link smearing methods currently in use, 
although an increased sensitivity to the smearing param- 
eter was found. 
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